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■ Abstract 

o- 

■ In this paper we consider the local and global well-posedness to the density-dependent incompress- 
ible viscoelastic fluids. We first study some linear models associated to the incompressible viscoelastic 
system. Then we approximate the system by a sequence of ordinary differential equations, by means 
of the Friedrichs method. Some uniform estimates for those solutions will be obtained. Using com- 
pactness arguments, we will get the local existence up to extracting a subsequence by means of 

■ Ascoli's lemma. With the help of small data conditions and hybird Besov spaces, we finally derive 
' the global existence. 

^' 



1 Introduction 



Elastic solids and viscous fluids are two extremes of material behavior. Viscoelastic fluids sho-w intermedi- 
ate behavior -with some remarkable phenomena due to their elastic nature. Their exhibit a combination of 
both fluid and solid characteristics and have received a great deal interest. It can also be regarded as the 
, consistence condition of the flo-w trajectories obtained from the velocity field u and also of those obtained 

■ from the deformation tensor F. Classically the motion of a fluid is described by a time-dependent family 
, of orientations preserving diffeomorphism X{t,x). Then deformation tensor F is defined as 

l/^ ■ 

• ^ dX(t,x) 

■ Applying the chain rule, we see that F{t, x) satisfies the following transport equation (see [T5]): 

dtF + U-WF ^Wu - F. 

. The viscoelastic fluid system of the incompressible in the Oldroyd-B model takes the following form: 

' dtU + u-Wu- fiAu + Wn = div{FF^), 
dtF + U- VF ^Vu - F, 
divu = 0, {u,F)\t=o^ [uii,Fo), 

where u is the fluid velocity, 11 is the pressure and F is the deformation tensor introduced above. Recently, 
the system has been studied extensively. Lin, Liu and Zhang in [T?], Lei, Liu and Zhou in [13_, Lin 
and Zhang in 16' proved the local well-posedness of (|l.ip in Hilbert space H'^ , and global well-posedness 
with small initial data. Local well-posedness can be proved by the standard energy method, while to 
obtain a global result, a very subtle energy estimate is applied to capture the damping mechanism on 
F — I. When one adds a linear damping term in the evolution equation of FF^ , which is the Cauchy- 
Green strain tensor, Chemin and Masmoudi [6] proved the existence of a local solution and a global small 
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solution in critical Besov spaces. We refer to [T^ and [19] for the well-posedness of the system in 
critical spaces. 

In fact, in the real world, the density usually depends on time. So we are more interested in the 
density dependent system, which is more close to the real situation. In this paper, we want to investigate 
the global well-posedness for the incompressible viscoelastic fluids described by the following system: 

at/9 + w • V/9 = 0, xeR^, t>0, 

dt{pu) + div(pM (g) u) - /iAw + Vn = div(det(F)-iFi^^), , , 

dtF + U-VF ^Vu - F, ^ ' 

divu = 0, {p,u,F)\t=o ^ {po,uq,Fo). 

The initial data {po,uo, Fq) are prescribed. 

Throughout this paper, we will use the notations of 

Ov 

(V^;),,, = (Vz;F),j = (V«),^fcFfe,,, (V • F), = d,F,,,, 

and the summation over repeated indices will always be understood. We also assume that oq = ^ — 1, 
Fq = Fq — I and Eq satisfy the following constrains: 

det(£;o + /) = !, dW{E^) = 0, (1.2) 

and 

dniEoij — djEQijn = EQijdiEQijn — EQirndlEQij. (1-3) 

Using these constrains, we obtain that 

det(£; + /) = 1, div(ST) -0, 
dmEij — djEijn — EijdiEijn — EijyidiEij, 

by Proposition 1 in [TS^. From the definition of F, we note that the assumption of det(ii^o + /) = 1 is 
nature. The first two of these expressions are just the consequences of the incompressibility condition 
and the last one can be understood as the consistency condition for changing variables between the 
Lagrangian and Eulerian coordinates. 

At this stage, we will use scaling considerations for (jl.ip to guess which spaces may be critical. We 
observe that ()1.1[) is invariant by the transformation 

{p{t, x), u{t, x),F{t, x), n(t, a;)) {pil^t, lx),lu{Pt, lx),F{ft, Ix), Pll{l^t, Ix)), 

{po{x),voix),Fa{x)) {poilx),lvo{lx), Fa{lx)). 

Definition 1.1. A function space E C (^'(R^) x J^'(R^))^ x (^'(R^))^><^ is called a critical space 
if the associated norm is invariant under the transformation (p{x) , u{x) , F (x)) — > {p{lx),lu{lx), F{lx)). 

Obviously ijf X (iJ^-i)^ X (iJ^)^^^ is a critical space for the initial data. The space H 2 however 
is not included in L°°, we cannot expect to get L°° control on the density and deformation tensor, when 
we choose Fq — I G (H^)'^^^ . Moreover, the product between functions does not extend continuously 
form X to , so that we will run into difficulties when estimating the nonlinear terms. 

Similar to the compressible Navier-Stokes system [7 , we could use homogeneous Besov spaces i3| i(IR.^) 

■ K 

(defined in [Tl, Chapter 2). B21 is an algebra embedded in L°° which allows us to control the density 
and deformation tensor form above without requiring more regularity on derivatives of po and Fq. Form 
now on, we define the density and usual strain tensor by the form 

1 



- 1,E:= F- I. 
P 



Then system (|l.ip can be rewritten as 



dta + u- Va^O, a; € E^, i > 0, 

dtUi + u ■ Vui - (a + l){pl\Ui - Vili) = Gi, , . 

dtE + u-VE^Vu-E + Vu, ^ ' 

divu = 0, {p,u,E)\t=o^{ao,UQ,EQ), 



2 



where = (a + l){djEikEjk + djE^j). 

Now we can state our main results. First define the following functional spaces: 

X^^C{[0,T];B^j')nL'{[0,T];B^y)xC{[0,T];Bl,), 

Y' =C{R+;B'^-^)nL\m+;B'2^^^) x C(M+; 5^'°°). 
Then the norm of and are defined by 

\\iu,E)\\x^ = ll"llz~(B=,-i) + ML^^(B';y) + 

\\{u,E)\\y^ = II"IIl~(r+,b=-1) + II"IIli(r+,b=_V) + II-^IIl-(r+,b;'°°)- 

Here B^ ^ denotes the nonhomogeneous Besov space and B^ ^ denotes the homogeneous space. The hybird 
Besov space B^^°° will be defined in the following section. 

Theorem 1.1 (Local well-posedness). Suppose that initial data satisfy the incompressible constrain 

N_ — — 1 — 

ao € B21, uq S B21 and Eq G B2i- Then there exist T > and a unique local solution for system 
ifTg)] with 

a e C([0,r];B2*i), {u, E) e Xj and VUe L^TiB^^i'^). 
Besides, the following estimate is valid 

\\a\\ N + \\{u,E)\\^ < C{\\ao\\ N + \\uo\\ + \\Eo\\ n), 
where C is a constant depending only on N and fi. 

Remark 1.1. We do not need the smallness condition on oq compared with the assumption of R. Danchin 
in ^ which consider the local well-posedness in homogeneous Besov space. The method was first intro- 
duced by R. Danchin in 1121 when dealing with the well-posedness of the barotropic viscous fluids in critical 
spaces. One can see in Section 3 that, for the technical reason, we could only study the local well-posedness 
on the nonhomogeneous Besov space without the smallness condition on oq . 

Theorem 1.2 (Global well-posedness). Suppose that initial data satisfy the incompressible constrains 
(r§ and ([13), ao G B^ , uq & B^^ and Eo G B^ ' with 

\\ao\\ + ||uo|| « _i + ll-Boll gf ,00 < A, 

where X is a small positive constant. Then there exists a unique global solution for system lll.2\) with 

aeC{R+;Bf'°°), {u,E)eY^ and VU e L\R+; B^^^^). 
Besides, the following estimate is valid 

<C(||ao|L«,^ + ||uo||.«-,+|li?o|| 

where C is a constant depending only on N and fi. 

The remained sections of this paper are structured as follows. In Section 2, we present some basic 
properties of Besov spaces. In Section 3, we will study some linear models associated to (jl.Sp . In Section 
4, the local theory for (II. 5p will be studied and the final section is devoted to discuss the global existence 
and to give the proof of Theorem 11.21 
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2 Littlewood-Paley decomposition results 

The proof of most of the results presented in the paper requires a dyadic decomposition of Fourier 
variables, which is called the Littlewood-Paley decomposition. The definition of Littlewood-Paley decom- 
position and Besov space were explained explicitly in [TJ [3J |3] . Here we state some classical properties 
for the Besov spaces. 

Proposition 2.1. The following properties hold true: 

1 ) Derivatives: we have 

C-'\\U\\^. < \\Vu\\^,-^<C\\u\\^^ . 

2) Algebraic property: for s > 0, B^ ^ n L°° is an algebra. 



3) Real interpolation: ^Bp^^, B^J^, 



A9si+(i-e)s2 



We recall from [T^ the following estimates for the product of two functions. Here we only give the 
results of homogeneous Besov spaces, and the same results also truth for nonhomogeneous Besov spaces. 

Proposition 2.2. The following estimates hold true: 

N 

\\uv\\.^^ N<\\u\\^sA\v\\^s2 if si,S2<— and Si + S2 > 0. 
2,1 2,1 z 

Let A = \/— A. For s e M, we denote K^'z — z), where z is the Fourier transform of z. 

The aim of this paper is to get the global existence of solutions to system (jl.Sp . For this, we define 
rf'-' = —K~^dju\ then = A~^djd^^ . Applying —A~^dj to the second equation of system (jl.Sp . we have 

dtd'^ - /lAd*^ + u • Vd - AEij = H. 

H will be determined in Section 5. Taking H as a function independent of d and E, combination with 
the third equation of system (jl.Sp , we have the following linear system 

dtE + Ad = R, 

dtd'-i - fiAd'^ - AEij =H. ^ ' 

Using the spectral analysis as in [7], we may expect that system p.ip has a parabolic smoothing effect 
on d and on the low frequencies of E, while expect a damping effect on the high frequencies of E. To get 
the optimal estimates, we need to introduce the hybird spaces which are defined differently for low and 
high frequencies. One can see the details in [7]. 

Definition 2.1. For /i > 0, r e [1, +00] and s S M, we denote 

\\u\\ 5.,. = J2 max{M, 2-^'-i\\K,u\\L^. 

Obviously we remark that ||it||j^s,oo k. p^^-i and ||it||j^s,2 — Also we need to introduce 

more accurate results which may be obtained by means of paradifferential calculus. It is introduced first 
by J. M. Bony in [2_. The paraproduct between / and g is defined by 

Tjg ^^Sq-xfi^qa- 
<?ez 

And define the remainder 

t9-p|<l 

We have the following so-called homogeneous Bony's decomposition: 

fg^Tsg + t.f + R{f,g). 

Now let us recall some estimates in hybird Besov spaces for the product of two functions which one 
can see Proposition 5.3 in [7]. 
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Proposition 2.3. Let r G [1, oo] and s,t gM.. There exists a constant C such 



that 



2 A'' A'' 

||T„^;|| < \\u\\ B-MUiy if s < min{l - - + y, y}, 

N 



\\R{u,v)\\ < M^Ms. if s + i>max{0,l-^}. 

^ (» 2,1 f 



Proposition 2.4. Let s,t gM.. There exists a constant C such that 

N 

< \\u\\gs^^\\v\\^ti, if s< — , 

N 

< ||u||5a,i||w||gt,oo, if s < — - 1, 

< if s + t>0. 

^2,1 

Proof. From the definition of T„t;, we can write 

KqfuV = Yl \{Sq'-luAg,v), 
\q-q'\<3 



whence 



\\Agfuv\\L2 < l|A,"u||i,<»||A,'v||i2 

\q'-g\<4 
q"<q'-2 

< 2«"*||A,.n|U2||Ag,t;|U2 

k;'-'7l<4 

q"<q'-2 

< Y max{/x,2-«"}||A5//M||i2 •max{/z,2-'''}-^||A5't>||i2 
x2«"^ max{/x,2-«"}-i • max{^, 2-«'}. 



|g'-9|<4 
q"<q'-2 



It is now clear that 

max{/i,2-«"}-imax{/x,2-«'} < max{l, 2«"-«'} < 1. 
So if s < ^, the convolution inequality yields 

gez 

For proving the second result, similarly, we notice that 

\\AqTuV\\L-< Y \\^q"U\\LA^q'V\\L- 

W-q\<i 

q"<q'-2 



< Y 2«"^||A,.«|U2||A,,^|U2 

\q'-q\<4 

q"<q'-2 

< Y max{At,2-«"}-^||A,»w||i2 •max{/z,2-«'}||A,,?;||i2 
X 2«"^ max{/x, 2-«"} • max{/i, 2-«'}-\ 



|9'-9|<4 
g"<9'-2 
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For q" < (7' — 2, we see that 

max{^,2-«"}max{/^,2"'?'}^i < 2'''-'?", 

then 



^2'^(^+*-f)||A,T„HU. <C^ J2 max{M,2-«"}-i2'"^||A,,m|U. 

9eZ |9'-q|<4 
q"<q'-2 

X max{^,2-«'}2«'*|lA^,w|li22(«-«')*2(«"-'?')(^-"-i). 



gez 9eZ|5'_q|<4 

q"<q'-2 



The convolution inequahty impUes 



gez 



2 

To prove the result on R{u,v), we note that 



ifs<f -1 



A,i?(u,'(;)= ^ (Aq^wA,--!;). 

This entails 

\\A,R{u,v)\\l2 <2'^f J2 \\Ag,u\\L2\\K,,v\\L^ 



q<q'+2 

<2«* max{A*,2-«'}|lA,-M|U2min{M-\2«'}||A,-i;|U2. 



9<9'+2 

If s + i > 0, then convolution inequality yields 



?6Z 



□ 



We will also use the so called Chemin-Lerner type spaces Lf^{Bp ^) which are described in detail in 
[5] . The case of nonhomogeneous Besov space can be defined in the same way. 

Definition 2.2. For p e [1, +00], s G K, and T e (0, +c»], we get 

and denote by ^) the subset of distributions u € S^'{0,T) x J?^^(R^) with finite \\u\\y^p ^ norm. 

When T = +00, the index T is omitted. We further denote Ct{BI^) = C([0, T]; S| n L^{Bl^) and 

3 The linearized equations 
3.1 The transport equation 

Here, we present a priori estimate for the linear transport equation which has been stated in [1] (Theorem 
3.14). 
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Proposition 3.1. Let 1 < p < pi < oo, r € [1, +oo], T > and ^ := 1 — -. Assume that 



s> ~Nmm{^,jr), if divu 7^ 0, 

s> -l-iVmin(^,i), */ divu = 0. 



Suppose ao G B^ ^., g G L^{0,T] B^ j.) and that a e L°°(0, T; solves 

dta + div(ita) = g, 
a\t=o = flo- 

Then there exists a constant C depending only on s,p,pi,r such that the following inequality holds, if 

ie[0,T], 



with 



V{t) = J^\\Vu{T)\\ « dr, if s<l + ^, 

V{t) = /o \\VuiT)\\s.-.dT, if s>l + ^ ors = l + ^,r^l 



-°2,1 



3.2 The linearized momentum equation 

When the density is close to a constant, we are led to study the following linearized momentum equations: 

dtu + v-Vu- fxbAu + 5Vn = /, 

divit = 0, (3.1) 
u\t=o = Wo, 

where 6 := a + 1 is bounded below by a positive constant b. That is inf b(x) > b. Before stating our 
result, let us introduce the following notation: 

At = 1 + 62^''"|1V6| 
for a e (0,1). 

Proposition 3.2. Let sG (1 — + '^^d, < a < 1. Also we assume a < if s > 1 and 

— " -1 

ao G B2i- Let uq be a divergence-free vector field with coefficients in B2 j. for r G [l,cx)], and f be a 

time- dependent vector field with coefficients in ^^^(iJj^^). m, u are two divergence- free time- dependent 

vector fields such that Vw G i^(0, T; B21) and u G C([0,r];B2/^ nL:^(B^+^). In addition, assume that 
^3.1\) is fulfilled for some distribution IT. Let Nq be a positive integer such that b^^ = 1 + SnoO satisfies 

inf 6jVo > 

Denoting fi := fi inf (a + 1), then there exists a constant C — C(s, N, pi, fi) such that if additionally, 

— xGM" — 

CA'^'^^Wa — S]y„a\\_ n < min{-fo, -r-uj, 
the following estimate holds for k = ^^j^ , 



< Ce^''^^^[\\uo\\s^^-.+A'i^{\\f\\^,^^^.-^^ 
with Vit)^J^{\\Vvir)\\. +22^'o||ap )dr. 



(3.2) 



^2,1 
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Proof. For positive integer Nq, we rewrite (|3.1I) as 

dtu + V ■ Vu — fibNo^u + 6Vn = / + i^AT^ , 

divu = 0, (3.3) 
u\t=o = uq, 

with i?7Vo = — SNoa)Au and bj\j„ — 1 + SN„a- 

Applying the operator to (13.31) . denoting Uq — AqU and 11^ — Agll, then we have 

dtUq + V ■ \/Uq — ^div(6jVoVUq) + VII^ 

= - Ag(aVn) + [v, Aq] ■Vu + Rq + A^Eno, 

where Rq := /iAg(5jV(,aAu) — /idiv(5ArpaVug). Muhiplying this equation by Uq, integrating by parts over 
yields 



W,\\h + tAl'^u.Wh < cWuqU^iWfqh^ + lib, A,] . Vu|U. 



^ (3.4) 

with fi :— fib. 

Denoting a := a — A_ia, since divu^ = 0, the last term of p.4p can be written as 

2/ Aq{a\7U)uqdx = 2 [ Aq{-Ty^n + Tv^a + R{a,\7n) + TA_,va^)uqdx 

by integration by parts and Bony's decomposition. According to Bernstein inequality, there exists a 
Ko > such that for all q > 0, we have Ko2'^\\Aqu\\]^2 < ||A^Vu||i2. Integrating over [0,T], p.4p implies 

\\Uq\\L2+fl2^'' f \\Uq\\L2dt 



<C\\AqUo\\L-+C J (\\fq\\L2 + \\AqEN„\\L-+tlS-l,q\\A-lu\\L^ 

+ \\[v,Aq] ■ Vu\\l2 + \\Rq\\L2 + 1 1 A, (Tysn) 1 1 ^2 + 1 1 A, (Ta _ , VaH) | U2 

+ ||A,(Tvna)||L2 + || A,i?(a, Vn)|U2)dt 
for all q > — 1- Elementary computations yield 

+ ||i?,|U2 + Ag]Vu|U2 + ||A,(Tvsn)|U2 + ||A,(Tvna)||L^ 



+||A,i?(a, vn)|U2 + WAqENoh^ + ||A,(rA_iVan)|U2j dtj 

7 

< C||uollB-i+Ail|w||£^(B. + i- = ) + ||/||£^(B.-i)+E/,. (3.5) 

1=1 

Now, we estimate the series of the right hand side of p.Sp term by term. 

As 1 — Y < ^ < ^ + '■^^^ Lemma B.2 in [TT] to bound 12- Indeed, 

i2 = {j2 2"^^'-'^ £\\[v,Aq]vu\\i2dty 



q>-l 

<C I \\Vv\\ .\\u\\J,. ^dt. 

Jo ^2.1 ^•'^ 



(3.6) 



8 



By Lemma B.3 in [TT], the following estimate 



q>-l 

< C/i|15jVoa|U^ «+„ I|m|Izi(b=+i-=) (^■'^) 

is proved. 

By the definition of Bony decomposition, we see that for any two functions / and g 
2«(^-i)||Tv/5||l^ <C 2«(^-i)||A,»V/|U22«"f ||A,,5||i. 

|g'-g|<4 
q"<q'-2 

<C l|A,"V/|U.2''"(f+"~i).2'''(^-i-")|lA,,Vg|U2 

|g'-g|<4 
q"<q'-2 

X 2('J~'j')(*~^)2^''"~*'''^"^~"^. 
Thus by the convolution inequality, /a is estimated by 

h < C||a- S'AToall « ||Vn||£i( . 1) 

1 / 

+ C||5Ar„a|| ||Vn||£,(^,-,-„), 

where we have used the above estimate with / — SnoO- and g = 11. 

Form Proposition 12.21 Theorem 2.82 and Theorem 2.85 in [1], /4 + /s is bounded by 

h + h < C\\a- SNoaW N \\S/n\\^^ ,g.-i^ 

As to /g, also the standard continuity result for para-product implies 

h <C\\a- SN„a\\^^^^fJMzi^(Bi+y (3.9) 

Obviously, Theorem 2.82 and Theorem 2.85 in T, implies 

h < C\\a- SNoaW N \\VU\\^, ,g. ^^ 

\^2,l) 

+ C||5Ar„a|| «+ II Vn||£^( . 

Thus, combining the above estimates for Ii to /g, we obtain 

< ||uo|Ib-^i + ||/|lzi(s.-i) +^^110- '5'A'oa|U^,„f Jl^llzi(s^^ 

\ 2.l) 

+ / (llVHI flkllBr/+A^2~""II^A^o«ll- Jl^llsn-O'^* (3.10) 
+ ||a-^Woa|l~ ,„f .l|Vn||£^(3.-i) +/fhllz^(s,=+^-=) 

V 2,1 / 

+ ll-^JVoalL^ «+„ ||Vn||£^(3. 
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On the other hand, VII solves the following elliptic equation: 

div(6jVoVn) =divL-f^jVo, 
with L = f + fmAu — v ■ Vu and Fjvq = div((a — SNga)\'lV). Apply to the above equation we get 

div(6jVoVn,) = divig - AgFNo + Rq, (3.11) 

with Rq = div(6jV(, Vllq) — Aqdiv(6jVo VII). Multiplying (I3.1ip by 11^ and integrating by parts, we obtain 

&||vn,||i. < (||divL,|U. + ||A,^^JU2 + |ii?,|U.)||n,|U.. (3.12) 

Bernstein inequality implies that 

fe2«||Vn,|U2 < c(\\divLq\\L2 + \\AgFN,\\L2 + \\Rq\\L2 + fo||A_in|| 



For — < (7 < , the second term can be estimated by 

9>-l 

Then from the assumption, 

C\\a~ SNoa\\ N < -b 
and Lemma B.l in [XT', for a G (0, 1), we have 

b\m\\Bi^ < IIQillBj.,, + \\SNoa\\ N^J^nw^.-. +6||vn||5.-o, (3.13) 

where 

Q = V(-A)"Miv, b^fi inf (SnoO- + 1) and - — < ct < — . 

x6R" 2 2 

If a satisfies that a < ct < by interpolation, we get that 

^l|vn||B,v < \\QL\\Bi,. + {\\SNoa\\ + 6)||vn||5.-. 

-°2,1 

< ||Qi||i3j,„ + (2^°"||5^„a|| « +6)||vn||g;i|vn||f.. 

-°2,1 

We conclude by Young's inequality and the estimate for the pressure in Proposition A.l in 11 that 

&l|Vn||Bj„ < A|||Si|lBj„. (3.15) 
Similarly, if a satisfies that — < f < — a, 

fe||vn||B.^, < \\ql\\bi^ + i\\SNoa\\ N +fe)||vn||^.+. 

^2,1 

y + o 



(3.14) 



<IIQi||Bj., + (||5^o«ll f +&)||vn||55j|vn||^,'. 

2,r 

Together with p.lSp . we can conclude that 

^l|Vn||Bj„ < A^f ||QL||bj„, 

if a < (7 < Y or < f'' < Therefore, we are led to estimate QL in L\,{B2~]^)- Since a < ^y^, this 
may be done by making use of Bony's decomposition, Lemma B.2 in ^llj and Propositions 2.2, 

+ /i||a - SNoa\\~^,_q.M^u\\i,^,gs-i^ (3.16) 
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Here we have used the fact that B^~j^ ^ B^'i'"" . Plus ((3?T6l) into (plO|) yields, 



< 



Ikolls-^i + ll/llzi (B-i) + Mt^^'IIo - ^A'oalU^^^f , Ikllzi (B=+i) 
+ f (llVHI « + m2^''"II«II „f 

+ / ||V«|| (3.17) 



The fact that 



^r(-S2i) 4 — 



and interpolation inequality imply that 



2,1 



Let X{t) — \\u\\j^^^^s-i^ + f^s+i-, + ||Vn||£i from the above estimate, we get 

rT ' , ' (3.18) 

+ / (||Vw|| N +22^«||a||°„ )x{t)dt. 



Then Gronwall lemma yields 

X(T)<Ce^^(^)(||uo|| 



Now we turn to give the priori estimate of the mixed linear system 

dtE + U-WE + Ad = F, 
dtd + u-Vd- fiAd ~ KE = G. 



(3.19) 



□ 



(3.20) 



We have the following proposition which one can see the details of proof in [7] . 

Proposition 3.3. LeA {E,d) be a solution of i3.20\) on [0,T] with initial data {Eo,do), 1 — < s < 1 + ^ 
and V{t) — ||Vu(t)|| . ivdr. Then the following estimate holds 



\\Eit)\\s.,^+\\d{t)\\^.-.+^, / (||£;(r)||5.,.+||d(r)||^...)dr<Ce^^W 

J 

X {\\Eo\\s...^ +\\do\\s^-^+fi re-^^M(||F(r)||5_ +||G(r)||^.-Orfr). 



4 Local well-posedness for data with critical regularity 

In this section, we will obtain the local existence of solutions to system p.Sp . We proceed by the following 
steps. 
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4.1 A priori estimates 

Let (ul, VII/^) solves the non-stationary Stokes system 



dtUL - fiAuL + VHl = 0, 
divMi, = 0, 

UL\t=0 = Uq- 



(4.1) 



It is easy to obtain that ul e C{[0,T];B^^^ ) n L^{0,T; B^^j^ ) and VHl G L^{Q,T\Bl^ 
that T has been chosen so small as to satisfy 



1-1 < ll'"o|| 



and 



where A will be determined later. 



Let M = M - itL, Vn = Vn - VHi, where (a, u, VH) satisfies ([T3)) on [0,r] 



Assume 



(4.2) 



(4.3) 



Suppose that 



a e C^{[0,T];B^\), u e Ci([0, T]; 53"^ n L^{B^y ) and VH e L\0,T;B^]^ ). We can deduce that 
(a, M, vn, £') satisfies the following system 



where 



dta + {u + Ul) ■ Va = 0, 

dtu + {u + Ul) ■ Vu - ^(a + l)Au +{a + l)Vn = F + G, 
dtE + u-VE ^ H, 
divu = 0, 
(a, u, E)\t=o = (ao, 0, i?o), 

F = — (u + ul)'S/ul + ^aAuL — aVIIi,, 

G = (a + i)(a,£;,fc£;jt + 9,i;,y), 

H = -ulVE + {Vu + Vul)E + Vu + Vul- 



(4.4) 



Denote Uq 



\\uo\\ Ul 
b7i 



inequalities are fulfilled for some suitable A, Uq and T: 



ivn. 



Assume that the following 



<2||ao| 



A^+^a-SNoa\ 



N < mm-^ — -0, „ riY, 



\E\\ 



L^{B^\) 



<mo\ 



(4.5) 



ivnii 



Li^iB^^ 



< XUo 



Then we are going to prove that they are actually satisfied with strict inequalities. Since (|4.5p depend 
continuously on the time variable and are satisfied with strict inequalities initially, a basic bootstrap 
argument insures that (14.51) are indeed satisfied for small T. For convenience, we denote 



U(T)^\\u\\_ N , N 



Ivnii 



'l\,{bZ ) 

First we prove (|4.5p 1 holds with strict inequality. Form Propositions 13.11 we easily obtain that 



L^iB^i) 



B^ 



^2,1 



(4.6) 
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where 

v{T)^ I (||vs|L« +||v,iLiL«)dt. 

If we choose A small enough such that 



-°2,1 -°2,1 



gC(-+^c/o) ^ 2, (4.7) 



then (I4.5l) i holds with strict inequality on [0,r). 
Similarly, form Propositions 13. Il we have 



\\E\\ . <e''^r.+i^''H\\Eo\\ . + f \\{Vu + Vu,){E + 1)\\ 



L¥iB2'i) ^ B^A Jo 

<6||i;o|| «, 

when Ai satisfies 

A A ~ 1 c( —+—Un) A A ~ 1 ,, ,, , . 

- + -t/o<3, e-^^-.-°)<2, - + -U,<-,\m,^. (4.8) 



Thus (|4.5|) ^ holds with strict inequality. 

According to the estimate (3.14) in T (page 134), we get that 

■ N 

II Aja|lL~(L^) < II Aj-aolk^ + Ccj2-i^ Ao\\V (u + ul)\\ n 

(-°2,l) 

where the norm of Cj equals to 1 and ^0 = 1 + llaoll " • By the definition of Besov norm, we see that 

||a-S'Ar„a||„ " ^ 2'?^||Aq(a-S'Ar„a)||ioo(i2) 



^ E E 2«^||AgA,a|Uo.(i2) 

j=No \q-]\<2 

00 







(-- 









2^"f ||A,ao||L^+CAo(- + -L^o 



Since oq G i?2^i and At < 2Ao, we can select Nq large enough such that 

C ^ 2^f ||A,ao||L2 < (2Ao)-('^+i) min{-5, — a*}. 

So (14.51) 9 holds with strict inequality provide 

CAo(- + -Uo) < (2Ao)-(^+i) min{l&, (4.9) 

Finally, we set T small enough such that 

4C22^°||aor„T<log2, 

B^^ 
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which combination with (|4.8p imphes 

eCV{T) < 4^ 

where V{T) is defined in Proposition 3.2. Thus from Proposition [321 we obtain that 

U{T) < Ce'^^(^)A^(|| - (u + ul)Vul +/iaAuL|l ^ «_i 

+ll«vni|| +||(a + i)a,i?,fei?,fc|| «_ 

+A*||a||- « I^lII +r||a||^ « \\E\\'i „ 



+r||a|L N {2\\E\\ N +1 +T||£;|L 
+ ||a|l jv llVnill N_, +2CAiA5,"+ir^||u|L 

^ ^,„^-iV (4.10) 



If we assume 



then we have 



A < ^M, (4.11) 



U(t) < 16cf(-[/o + 2-||ao|| «) + 2||ao|| n)X + 36T\\Eq\\^ j, 

V /i fl B^, B^, B^ 

+ 144T||ao|| JV ||i;of « +48r||ao|| N\\Ea\\ jv+6r||i;o|| m 



3 "2" 

2,1 



+2r||ao|| «) +64C2^°"/^-(l + ||ao|| nT°'+^T'^\Uo 

B^/ /I B^ 

< CoiUo + l)X + CoT + Co2^'"'TXUo, (4.12) 

where Co is a general constant depending only on ||ao|| «, ||ao|| jv , /i, fi. Hence, selecting Uq = 

8Co([/o + 1), for fixed A which determined by (14. 7p . (j4.8p . (|4.9I) and (j4.1ip . we can choose T smaU enough 
such that 

Co2^°"r<i 4C22^°|laoPivT<log2, CqT < ^AC/o. (4.13) 

4 R 2, 8 



This implies (|4.5P a holds with strict inequality. 

4.2 Priedrichs Approximation and uniform estimates 

Let be the set of functions spectrally supported in the annulus C„ = G |^| < n}. Jn denotes 
the Friedrichs projector maps to L^, defined by 

^J,,u{S) = 1c„^m(C) for aU i e M^. 
We aim to solve the system of ordinary differential equations 

= Fn{a,u,E), 
iu = G„ia,u,E), 

iE^H„{a,u,E), ^ ■ > 

{a,u,E)\t=o = (J„ao, 0, J„i?o), 
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in Ll X [Lir X {Lir with 

Fn{a,u,E) = - Jn{u ■ Va), 

G„(a, u, E) — — Jniu ■ Vu) + /iJ„ (6A?2) — J„(aVni) 

+ fiJniaAuL) + Jnibdiv{EE'^)) + JnibdivE) 

+ J„(6'Hb(- J„(w • Vii) + fiJ„{aA{u + ul))) 

+ JJbHbiaVUL)) + JnibHb{bdiv{EE^) + bdivE)), 

Hnia, u, E) = -Jn{u ■ VE) + J„(Vw • E) + JnVu. 

Here u = u + ul, ul is the solution of (I4.ip . Hb denotes the Hnear operator F i— VII, i.e. VII = Hb{F) 
is the solution of the elliptic equation 

div(feVn) = divF. 

The map 

{a,u,E) I — > {Fn{a,u,E),Gn{a,u,E),H„{a,u,E)) 

is locally Lipschitz with respect to the variables {a,u,E). Then we can conclude that the ordinary 
differential equations has a unique solution (a", m", i?") in the space C^([0, T*); L^). T* is the maximum 
existence time of {a^,u",E"-). Then using the elliptic equation we can get the existence of VII" £ 
CH[0,T*);Ll). 

Now we want to prove that T* may be bounded from below by the supremum T of all the times satis- 

— 

fying (|4.13l) . and that (a", m", E") is uniformly bounded in Xr^ . Since J„ is an orthogonal projector, 
it has no effect on the priori estimates which were obtained in Section 4.1. Hence, the priori estimates 
applies to our approximate solution (a", u", E", VH") which independent of n. And the estimate in (|4.5p 
to (a", w", i?", VH") ensue that it is bounded in i°°(0, T; i^J. So the standard continuation criterion 
for ordinary differential equations implies that T* is greater than any time T satisfying (j4.13p and for all 
n > 1, 

||a"|L N < 2||aoll N , 

A-||a"_.,,«"||^^^^^^„^^<min{±5, ^,}, 

||m"IU iv+1 +||vn"||„ < ac/q. 

4.3 Compactness arguments 

We now have to prove the convergence of (a", w", i?"). This is of course a trifle more difficult and requires 
compactness results. Let us first state the following lemma. 

Lemma 4.1. (a" ,u" , E") is uniformly bounded in 

Cn[0,r],i?2^i ')xC^([0,T],B2^i ' + ^2^ ')xCH[0,TlBl, 

for N>3. 

— —1 

Proof. We first prove that dta" is uniformly bounded in L^{B2i ), which yields the desired result for 
a". 

We observe that a" satisfies 

9ta" = -J„(m"- Va"). 

According to the uniformly estimates in Section 4.2, m" is uniformly bounded in L'^(B2i) and a" is 
uniformly bounded in L^{B2i)- The definition of u£ obviously provides us with uniform bounds for it 
in L|n(i?2^i)- So we can conclude that 

||9ta"|| N , <cf||u"|| N +IKII « )||a"|| 



(4.15) 
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— — 1 

which imphes that dta" is uniformly bounded in L'^{B2i )■ 

——1 

Similarly, we show that dtE" is uniformly bounded in L^{B2i )■ Let us recall that 
dtE'' = J„ ( - (S" + ul) ■ VE'' + (Vu" + \7ul) ■ E'' + (Vm" + Vu^) 
By the continuity of Paraproduct in Besov spaces which is stated in Proposition 12. 1[ we have 

\\dtE"\\ <cf||u" + <|| «)f||S"|| ^+1), 

——1 

which implies that c?ti?" is uniformly bounded in Ll^^B^i )■ 

— —2 — — 1 

Now we turn to prove dfU^ is uniformly bounded in L'^{B2i ) + L'^{B2i )■ Note that dtu"^ satisfies 



dtu" = Jn[Gn - Vn"(l + a") + A*Au"(l + a") - VHJ 

- (w" + O • Vu" + A'a"A< - + • V<) , 
where G„^i = (a" + '^){djE^^.E^^. + djElj). As the above estimates, we know that 

^Au"(l + a") - (u" + ■ Vu" + ^a"Aw2 - + O ■ "^u'l + VHJ 

— —2 

is uniformly bounded in L^{B2i )■ By the expression of Gn,i, we have 
||G„|| N_, < CrUlla" 



which implies that G„ is uniformly bounded in L^{B2i )• Now we devote to estimate VII". We split 
Vn" into Vn^' and VB^, and their satisfy 

div(6"VB5') = divG„, div(6"Vn'2') = divF„, 

where En ~ ^Au"a^^ — (m" + u£) • Vu" + /Lta"Aw2 — (w" + u£) ■ Vu£. By the estimate for the elliptic 
equation in Proposition 13.21 we get 



||vn"i| N , <c\\Gj « i,||vn"i| « , <g||f„|| « 

From the above discuss, we know that \\Gn\\ « , and IlKJI « „ are bounded, and VB" is 

N__2 ——1 — — 2 — — 1 

uniformly bounded in L|.(i?2^i ) + L'^{B2i )■ So (?tu" is uniformly bounded in L|n(i32^i )+L^{B2i )■ 
Thus we have proved the lemma. □ 

We can now turn to prove the existence of a solution. The procedure is similar as been used in [7]. 
We can see in [7] that the approximation solutions are convergence in the term of subsequence by Ascoli's 
theorem. So we omit the details. The same argument to iV = 2, we can also prove that dtu" is uniformly 
bounded in L^[0,T], B2 I). 

4.4 Uniqueness 

Assume that we have two solutions of (jl.Sp . (ai, ui, E'l, VBi) and (02, U2, E2, VB2) with the same initial 
data satisfying the regularity assumptions of Theorem ll.il We first consider the case A'' > 3. 
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Set ai-a2 = Sa, m - U2 = Su, VIIi - = VSU, Ei - E2 = 5E. Then {5a, 6u, V5P, 5E) satisfies 
tlie following system 

dfSa + U2 ■ WSa = —Su ■ Vai, 

dtSu + U2 ■ VSu - ^(1 + ai){A6u - VJH) = 5G + SH, 

dtSE + U2-VSE = SL, (4.16) 
div 6u = 0, 
[ {Sa,Su,6E)\t=o = (0,0,0), 

where 

6H = —Su ■ Vui + fj,SaAu2 — SaVIl2, 
6Gi = (ai + l)dj5EikEi,jk + (oi + l)djE2,ikSEjk 

+ SadjE2,ikE2^jk + 6adjE2Aj + {a-i + l)djdEij, 
SL = -Su ■ VEi + Vu2 ■ 5E + V6u ■ Ei + V6u. 
From Proposition 13. 1[ we have 

+CT\\SNoaiL^^N^„Ti\\Su\\^ \\6u\\''^ „ , 

where V{T) — ||Vm2|| n and 77 e (0, 1). Similarly, we have 



Next, denoting = II (SiilL « , « +||V(5n|| n „, the estimate for linear momentum 

equations. Proposition [32] guides us to get 

SU < e'^^C^^^^yf / \\Su\\ N_^\\\/ui\\ N +fi\\Sa\\ iv_i||Au2|| n _^dt 

' ^Jo B^^ B^^ B^ B^ 

+ r\\6a\\ ._,\\VU2\\ .^,dt+\\SG\\ . 

JQ B21 -D2.1 -'^r(-°2,i I 



+Mi,t||<5u|| «_ , (4.19) 
with A\^T — 1 + 62^°"||ai|| « and a G (0, 1). By interpolation, the last term can be bounded by 



' b'^ 

-°2,1 



Young's inequality and the uniform bounded of the solution imply that 

su<Ct f {i + \\u2\\ «+i + ||vn2|| «+j 

Jo B^i -B2.1 

{\\5u\\ ._2 + \\5a\\ ._,dt+\\5G\\ « ) 



dt. 
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By the expression of 5G, we have the following estimate 

\\5G\\ 



< I (||<5a|! «_j|ii;2||% +(||i?i|| « +||i?2|| 

Jo ^ -^2.1 -BsTl -^sTl ^2,1 ^2,1 

+ IIHL f-.p2|| « «||<5i?|| 

2.1 ^2.1 ^2.1 ^2,1 

+ ||ai|| « +P2II «)||<5i?|| 

^2,1 ^2,1 ^2,1 ^2,1 



Combination the above estimates, we know that 

||<5a|L jv_, + UE\\_^ +5U 



<Ct f (1 + IIU2II «^, + ||vn2|| 

Jo ^2,1 ^2,1 

X (ll^ull + ||(5a|| + N_Adt. 



which yields 

||(5a|L « 1 = y£;|L iv 1 = (5C/ = 0, 

for small enough T. A standard continuity argument allows us to know the uniqueness on [0,T*), T* is 
the lifespan of the local solution. This finish the proof of the uniqueness of Theorem 11.11 when > 3. 

In the case of A = 2, the above proof fails because — 1 = 0. Hence we may be tempted to estimate 
{Sa,6u,SE,VdU) in 

Now we give the details of the proof. From Proposition 13.11 we have 

/ )\ (4-20) 

When T is finite, ||(5w||£i ) < W{T), where W{T) is finite. Hence 

/ WiT) \ 

-j- (4.21) 

-^'r(-°2,oo) 

Also we can get 

^^^^h^iBi^) 



< Ct{\\5u\\^,^^^,^^+ \\u2\\BiJSE\\soJt) (4.22) 
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W{T) 

Then Gronwall inequality implies 

W{T) 



dt. 



/ W(T) \ 

W^^'hnBi^)) ^ ^^W^^huBi^) i°g + ||fa|u,^ . ) • (4.23) 



Next, denoting SU = \\Su\\j^ao f^^-i ^ +^i\\Su\\y^i ^ + \\WSII\\y^i ^g-i ^ the estimate for linear momentum 
equations, Proposition 13. 21 guides us to get 



||(5u||^-i WWuiWb^ + ^J.\\Sa\\Bn \\u2\\B^dt 

+ ||<5a|Uo^||Vn2bo/t+||5G||£.(^-^) (4.24) 

+ ^iAi^T\\H\LUBlT))^ 

with Ai T = 1 + ^2^°"||ai||roo^Ri -1 and a E (0, 1). By the expression of SG, with the same calculus as 

' T ^ 2 ,1 ) 

N > 3, we have the following estimate 

||'5G||£^(5-^) < Ct f {\\5a\\Bi^ + \\6E\\Bi^)dt. (4.25) 
^ 

Hence combination the above estimates, together with interpolation and Young's inequality, we know 
that 

5U{T) < CtJ^ ((1 + \\u,mB2jSUit)\og (e+ ■ffi))^^ (4.26) 

which yields dU = on [0,T] by Osgood Lemma (Lemma 3.4, Chapter 3, jl]) because 

/ ^dr = +00. 

•^0 r log i^e+^j 

This finish the proof of the uniqueness of Theorem 11.11 when N = 2. 



5 The Global Theory For Small Initial Velocity 

In the above section, we have proved that there exists a unique local solution (a, u, E) of (II. 5p in 

N_ ^ — ^ 

C{[Q,T];Bl-^) X . We have used the estimate for VH in ((3T3l) . That is the reason why 

we work on the nonhomogeneous Besov space. We rewrite p.l3p as follows: 

^livniisj, < ||ql||b|, +||a|| « livniis, 

B. 



2ir-oll~ "-"2,1 " "nT" II -"2,1 

2,1 



While fortunately, the assumption on ao in Theorem ll.2l can avoid the estimate of VII. More precisely, 
the second term can be absorbed by the left hand side due to the smallness condition on a. Thus using 
the same method as in Theorem 1 1.11 we obtain that there exists a unique local solution (a, u, E) of p.Sp 
in 

C{[0,T*y,B,\) X L'{0,T*;B,-/')nC{[0,T*);Bl, x C([0, T*); B^^J, 

where T* is the maximum existence time of {a,u,E). From the assumptions in Theorem 1.4, using 
Proposition 3.1, we can easily obtain 

. jV _-! . N__'i 

aeC{[0,T*)-Bl^ EeC{[0,T*)-Bi, '). 
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We define d*-' = — A ^djU^, then = A ^djd^^. Applying —A ^dj to the second equation of system 
(Hm), we get 



+ A-^dj(u ■ Vu' + (a + l)d.,U - ^laAu' - G J , 



(5.1) 



where Gi ^ {a + l){djEikEjk + djEij) is defined in Section 1. Note that the compatibihty condition 
p.3p . we have 



= A~^dk{dkEij + EikdiEij - EijdiEik) 
= —AEij + A^^dk{EikdiEij — EijdiEik). 

Combination ([ET|) . with ([T3|) yields 

dta + u ■ Va = 0, 

dtd'^ + u ■ Vd'^ - fiAd'^ - AE,j = H, 
dtE + U-VE + Ad'3 = R, 
d'^ = -A-^dju\ 
divu — 0, 

(a,M, £;)|t=o = (ao,uo,-Bo), 



(5.2) 



(5.3) 



where 



H = -u- V{A~^dju') + A^^dj [u ■ Vu' + (a + l)diYl - fiaAu' 
- (a + l)diE,kEik - adiEa) - A-^dk{EikdiE,, - Ei,diE,k), 



Denote a = ||ao|| + |luo|| . « _ 



i? = dku'Ekj. 

I-E'oll .oo • We are going to prove the existence of a positive M 

]3 



such that, if a is small enough, the following bound holds 

l!a|U -iv + \\iu,E)\\ N < Ma. 

This estimate is the direct product of the following proposition. 
Proposition 5.1. // 

||a|L + ||(u,i?)|| " < 2A/a, T G (0,r*) 



(5.4) 



Lf(B^^ ■ ) 



Y 2 



then, we have 

\\a[ ^ 
when a is small enough. 

Proof. First, from Proposition 3.1 in f5^, we obtain 



a\\^^^^N^^^ + \\{u,E)\\ N <Ma, 



L^iBZ) ~ 



2,1 



where V{T) = ||Vw|| . « dt. If we assume a small enough such that 



2,1 
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then we have 



for Ad = 4Ci. From Proposition 3.4, we have 

YrJ' \ B,r BJ. 



R ^ 

(5.5) 



+ ||i?|| +||i?|| 

We want to bound ||i?|| -.n ^ and \\H\\ . w_i • With the help of Proposition 12. 3[ we have 

\\R\\ <C\\E\\ ^iv^||Vw|| jv < CAf^a^ 

Now, we devote to estimate \\H\\ n . From the expression of H, the trouble is the estimate for 
Vn. Applying div to the momentum equation of (|1.5p yields 

d,{{a + 1)8,11) = di{-u ■ Vw* + naAu' + L,), 

with 

= (a + VjdjEikEjk + adjEij. 
Here we have used div(i?^) = 0. Then by the estimate of elliptic equation, the following bound holds 

||Vn|| «_! < CllulP JV +C||a|| N \\u\\ + \\Li\\ iv_i . 

Note the expression of Li, we only need to estimate ||a9j£'ij|| ^ . n_i . In fact, by Proposition 12 .41 with 



s — — Y ~ have the following 



||aa,-£;,,|| N , < C\\a\\ ^ li-Bll ^ 



Others in Li are estimated similarly. Hence we get 



(5.6) 



IVHII < CM^a^(l + Ma). 

l't(bZ ) ~ 



For the term of H, using the estimate of VH and Proposition 12. 2[ we have 



li^ll <C\\uf .„ +CM||a|L , ~ l|"ll 



+ C||a|| « ||VH|| +C||VH|| 

"L^iB^y L}^{bZ ) iT(S2"^i ) (5.7) 

+ C||(a+ 1)£;- V£;|| iv 1 + \\aVE\\ 

L},{B,y') L^riBi^ ) 

< CM^a^{l + Ma). 

Plugging the estimates on H and R into (|5.5p . noting that ~ —A^^dju^ , we have 

llalL ^ + ||(u,i;)|| « < Cae^A'^" fa + CaM^a^ + CaM^a^) 

iT(s.^' ) V " ^ ^ (5.8) 

< Ma, 

when Af = 46*2 and a satisfies 

g2Ma < 2, 26*2^ Af a < i, 2CH'I^a^ < -. 
~ ^ ~ 4 ~ 4 

Then, we finish the proof of Proposition 15 . II for M = max{4Ci, 4C2}. □ 

Now we can give the proof of the global existence. From the standard continuation method and 
Proposition 15. 1[ we easily obtain that (|5.4p holds. Combining the local existence, if T* is finite, then the 
lifespan of the solution is greater than T* . Hence T* ^ 00 and we finish the proof of Theorem 11.21 
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